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Spin Damping Monopole 

Akihito Takeuchi* and Gen Tatara 
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We present theoretical evidence that a magnetic monopole emerges in dynamic magnetic 
systems in the presence of the spin-orbit interaction. The monopole field is expressed in 
terms of spin damping associated with magnetization dynamics. We demonstrate that the 
observation of this spin damping monopole is accomplished electrically using Ampere's law 
for monopole current. Our discovery suggests the integration of monopoles into electronics, 
namely, monopolotronics. 
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The magnetic monopole predicted by Dirac in 1931 1 is a unique particle that arises from 
singularity. 2 In high-energy physics, a monopole emerges if the electromagnetic interaction 
in the world, described by a U(l) algebra, is a result of the symmetry breaking of a unified 
force having a higher symmetry of SU(5). 3,4 Intensive effort has been exerted to find a grand 
unified theory (GUT) monopole by waiting for a monopole created in the early universe to 
go through superconducting detectors 5 and by detecting its ionization; 6 however, no evidence 
has been found thus far. The energy needed to create a GUT monopole is about 10 17 GeV, 
and so creating one in an accelerator on earth is impossible. 

Since a monopole is a consequence of symmetry breaking, it exists in solids too. Symmetry 
breaking in solids occurs at an energy much lower than the GUT energy, typically below 1 eV, 
and therefore experiments on this are feasible. The most well-known monopole in solids is the 
hedgehog monopole (HHM), which arises in magnetic materials; 7 it is the SU(2) counterpart 
of the GUT monopole. The key interaction for the HHM is the coupling of the conduction 
electron to the magnetization represented by the vector M(r,t), which depends on the space 
coordinate r and the time t. The electronic spin, represented by a vector <r, is polarized by 
M owing to the coupling 

H sd = -JM ■ <x, (1) 

where J is the coupling constant. Since M is an external field for electrons, the SU(2) symme- 
try of the electronic spin is broken. By diagonalizing H s & and choosing the spin quantization 
direction to be along the z-axis, the electrons are described as spin-polarized and interacting 
with a SU(2) gauge field, 7 (// = t, x, y, z and a = x,y,z are the indices in the coordi- 
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nate space and spin space, respectively). In the adiabatic limit, i.e., when J is large, only the 
z component of the gauge field, A*, survives and it acts as a U(l) gauge field. The HHM 
originates from the deviation from the perfect U(l) symmetry, i.e., from the perpendicular 
components A* and A^. In fact, the electrons feel the SU(2) gauge fields whose strength is 
F^ v = d^A® — d v A°^ — fie/fyeabcA^A^, where e is the electric charge, h is the Planck constant 
divided by 2ir, and e a b c represents the antisymmetric tensor. When this field is projected onto 
the U(l) space, we obtain the electromagnetic field tensor as F* v = d^Af, — d v A z + 
where = —(2e/h)(A^Af / — A^Af,). The anomalous field strength ^j,, which in terms 
of n = M/\M\ reads <3? Ml , = —(h/2e)n ■ (d^n x d u n), represents the HHM. In fact, the 
field strength satisfies (l/2)e flvap d l ,F* = 0, whose components read V x E + B = —jh, and 
V-.B = ph) where the electric and magnetic fields are Ei = —V{A\—A\ = — (h/2e)n-(hx V^n) 
and Bi = eij^jAt = (h/Ae)e,ijj t n ■ (Vj-n x V^n), respectively. The monopole current (jh) 
and its density (ph) for the HHM are given as j^i = —(3h/4e)€iji ; n ■ (V^n x V^n) and 
Ph = (^-/4e)ejjfcVjn • (V^n x V^n), respectively. Although the HHM is mathematically al- 
lowed, experimental realization has not been achieved thus far. In fact, as seen from the 
expression for jh and p^, the HHM disappears in common magnets where the length of the 
local magnetization, n, is constant. In addition, the boundary condition at infinity for the 
HHM would not be easy to realize experimentally. 

In this paper, we search for a different monopole in magnets, which exists in conventional 
ferromagnets where the local magnetization length is constant. Such a monopole current 
creates the rotational electric field via Ampere's law. This means that the monopole current is 
an anomalous angular momentum source that induces the rotational motion of electric charges. 
To realize such a monopole in magnetic systems, a coupling between spin and electron orbital 
motion is, therefore, essential. Such a coupling is known to emerge from a relativistic effect, 
namely, the spin-orbit interaction. The spin-orbit interaction exists in all elements including 
magnetic ones and is particularly strong in heavy elements such as platinum and gold. 

Our aim in this study is to prove the existence of the above monopole theoretically. 
Since the spin-orbit interaction explicitly breaks the SU(2) invariance, we cannot follow the 
derivation of the HHM shown above. Instead, we will directly calculate the effective electric and 
magnetic fields for the electron spin based on a nonequilibrium Green's function formalism, 
and derive Maxwell's equations they satisfy. 

We consider two types of spin-orbit interaction. The first is that from a uniform field, £Jr, 
namely the Rashba interaction. 8 Such a field is realized at interfaces and surfaces. 9 The second 
is that from a random potential, v\, induced by heavy impurities. The spin-orbit interaction 
thus reads 

H so = -^{X R E R -X i Vv i )-(pxa), (2) 
where p is the electron momentum and A is the coupling constant (the subscripts R and i 
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characterize Rashba and impurity- induced spin-orbit intereactions, respectively). The inter- 
action with the magnetization is described by -£f s d- The Hamiltonian of the present system is, 
therefore, given as H = (p 2 /2m) + V{ + H S( i + H so , where m is the electron mass. The magneti- 
zation we consider in H s d is dynamic. Dynamic magnetization, when coupled to the spin-orbit 
interaction, generates an electric charge flow. 10-12 The pumped electric current j is calculated 
by evaluating a quantum field theoretical expectation value of the electron velocity operator 
v = — (z?i/m)V + (l/h)(X R E R — AiV^i) x a. By using field operators for electrons, c 1 and c, 
the electric current thus reads j = —etr(c'vc), where tr denotes the trace over spin indices and 
the bracket represents the expectation value. It is written in terms of the lesser component of 
the nonequilibrium Green's function, 13 defined as Gf s ,(r,t;r',t r ) = (i/h)(c s ,(r f ,t')c s (r,t)) (s 
and s' are spin indices), as 

ti 2 

j(r,t) =etr({ — (V r - V r ,) + i[\ R E R 

- XiVviir)} x *}G< (r, t; r>, t)) r , =r . (3) 

This quantum field theoretical expectation value is evaluated by solving the Dyson's equation 
for the nonequilibrium Green's function defined on the Keldysh contour (C), 

G ss >(r, t; r', t') =5 s ^g s (r, t; r' , t') 

+ J d\" J dt"g a (r,t;r",t") 

x JM(r",t")-* ss „ 

+ i{[\ R E R - \iVvi(r")] x V r »} • <r ss ») 

xG sl/s/ (r"X;r',t'), (4) 

where G ss /(r,t;r',t') = —(i/h){Tc[c s (r,t)cl,(r',t')]) (Tq is the path-ordering operator) and 
g denotes the free Green's function. 

In the calculation, the impurities are approximated as random point scatterers, and aver- 
aging is carried out as (vi(r)vi(r'))i = niuf5 3 (r — r r ) (ni and u\ are the impurity concentration 
and the strength of scattering, respectively). 14 The impurities give rise to an elastic lifetime 
for the electron, r, which is calculated in metals as r = h/2imiU 2 v (y is the density of states 
per volume). 15 The Dyson's equation is solved by treating A and J perturbatively to the first 
and second orders, respectively. We consider a sufficiently slow dynamics of magnetization, 
namely, Qt <C 1 (f2 is the frequency of magnetization dynamics), and assume that the mag- 
netization structure varies smoothly in the space compared with the electron mean free path 
£, i.e., q£ <C 1 (q is the wave number of magnetization profile). The leading contribution in 
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this case turns out to be 16, 17 

k,k',qi,q2 U!,Cli,Cl2 

x fj(M qi , ni x M q2 , Q2 ) x [^fE R \glJ 2 

Aft 2 A; r o i r i4i 

+ j^ili + <l2)£k\g k ,J \g k >. 



,0J\ 



-DVp(r,t), (5) 

where V is the system volume, D = 2£ft/3?ti is the electron diffusion constant (sf represents 
the Fermi energy), f u is the Fermi distribution function, g T is the retarded Green's function 
defined as g T k = [ftw — e k + (ih/2r)]~ 1 , and = ft?k 2 /2m. The last term is the diffusive 
contribution arising from vertex corrections, where the electric charge density p is 19 



X 



h 2 v 

e -iq-(r-r')+in(t-t') 



EE 



Dq 2 T + iVtr 



x{E R x[M(r',t')xM(r',t')]}. (6) 

Summing over the wave vectors and frequencies in eq. (5), the electric current is obtained as 

16e^Ai«/ 2 £FT 2 , 



3ft 2 
4eis\ R J 2 T 2 



-V x (M x M) 



E R x(M x M) - DVp. (7) 



^ 2 

This result is rewritten using the effective electric and magnetic fields, E s and B s , respectively, 
as (p and a c are the magnetic permeability and electric conductivity, respectively) 

j = -V x B s + a c E s - DVp, (8) 
A* 

where the effective fields are defined as 20 ' 21 

E s = -a R E R x N, 

B s = -fcN. (9) 

Here, N = M x M is a vector representing the spin damping torque (inset in Fig. I). 22 The 
coefficients a R and /?; are a R = 4ev\ R J 2 T 2 /a c ft 2 and /3; = 16evp\i J 2 EFT 2 /3ft 2 , respectively. 
The effective fields calculated here are those acting on the electronic spin in the same manner 
as the effective fields from the HHM. Clearly, the fields [eq. (9)] do not satisfy Faraday's law 
or Gauss's law of conventional electromagnetism, but those with monopole contribution: 

V xE s + B s = -j m , 

V B s = p m , (10) 



4/9 



J. Phys. Soc. Jpn. 



Letter 




Fig. 1. (Color online). Schematic illustration of monopole pumping and detection in thin ferromag- 
netic film attached to a nonmagnetic layer. Magnetization (M) precession is induced by applying 
an oscillating magnetic field. The Rashba field exists at the interface and creates the monopole 
current j m near the interface. The width of the monopole current distribution, d, is comparable to 
the decay length of the magnetization at the interface. The monopole current induces the electric 
current j via Ampere's law at the interface. The impurity spin-orbit interaction directly induces 
positive (+) and negative (— ) monopole charge distributions, p m , at the two edges. This monopole 
distribution generates an electric current again perpendicular to the monopole current, as is seen 
from eqs. (8) - (10). Inset: depiction of the spin damping vector N = M x M arising from the 
precession of magnetization. The component of the spin damping vector perpendicular to the 
precession axis vanishes when time-averaged, leaving N along the axis as the DC component. 



where the monopole current and monopole density, respectively read 

jm = "rV x (E R x N) + ftTV, (11) 

and 

Pm = "AV • N. (12) 

Equations (10) - (12) are the central results of this paper. We have thus proved that 
a monopole exists when spin damping occurs, namely, a spin damping monopole. The spin 
damping monopole is a composite object made from a magnetization configuration in the 
same manner as the HHM. It satisfies the conservation law p m + V • j m = 0. 

Our result obtained in a disordered system is a general one and can be extended to a 
clean system. In fact, the same monopole exists in a clean case, but only the coefficients cur 
and fa appearing in the monopole density and current are changed. The same applies to the 
HHM; when we take into account the third-order contribution in J, our analysis correctly 
reproduces the HHM, which was discussed in a clean limit only. Our approach is, therefore, a 
novel method of identifying monopoles. 

The spin damping monopole is unique since it does not require a particular non-coplanar 
spin structure like a hedgehog, and so it exists quite generally in magnetic systems. The 
simplest candidate for creating the monopole would be a thin ferromagnetic film put on a 
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nonmagnetic material, as shown in Fig. 1. We choose the z-axis perpendicular to the film. 
A Rashba-type spin-orbit field would then arise at the interface along the z-direction. 9 We 
excite the precession of the uniform magnetization by applying the alternating magnetic 
field in the yz-plane in the presence of a static field along the x-axis (ferromagnetic reso- 
nance 22 ). The precession results in a spin damping vector with a finite time average, N, 
along x-direction (inset in Fig. 1). In the present case with uniform magnetization, spatial 
derivatives in eqs. (11) and (12) arise at the interface and edges, where the magnetization 
vanishes. The Rashba interaction contributes to the DC monopole current at the interface 
as Jmi = — a RE-£i(dN /dz) ~ —(an/d)EnN, where d is the spatial scale of the magnetiza- 
tion decay at the interface. The monopole current driven by random spin-orbit impurities, 
on the other hand, vanishes when time-averaged. The total DC monopole current thus reads 
j m = —e x (aji/d)E^iN (e x represents the unit vector along the x-direction). This monopole 
current at the interface generates an electromotive force along the y-direction via Ampere's 
law for the monopole. The monopole density induced by the random spin-orbit interaction 
arises at the edge of the ferromagnetic film since V • N ~ dN x /dx is finite there. The induced 
monopole density at the two edges is = ^f((3i/d)N , where the sign is positive on one side 
of the edge and negative on the other side. The monopole then produces a magnetic field 
along the x-direction as B s = —e x f3{N. This field creates an electric current in the y-direction 
via the conventional Ampere's law. The total electric current density generated by the spin 
damping monopole [eq. (8)] thus reduces to j = — e y [cr c anEn + (f3i/[id)]N. 

When spin damping arises from the magnetization precession with a frequency 
f2 and an angle 9, the monopole-induced current density is estimated as |j| = 
(eklQ/7r 2 )(jT/h) 2 sin6»[(A R /e F ) + (A/3k F d)(A i /v 1 )], where A R and A ; are the energy of the 
Rashba and impurity spin-orbit couplings, respectively (hp is the Fermi wave vector). In dis- 
ordered ferromagnets, J/ep ~ 0.1, epr/h ~ 10 and k^ 1 ~ 2 A. The Rashba interaction 
can be enhanced on surfaces and at interfaces, resulting in A r /£f ~ 0.1 (Aj/ep is gener- 
ally smaller). 23 When 6 = 30° and Q = 1 GHz, the electric current density is thus 2 x 10 7 
A/m 2 which is sufficiently large for experimental detection. In addition to DC, there is an AC 
component in eq. (8), which would be accessible by time-resolved measurement. 

We note that the electric current estimated here is an initial current that arises when 
the pumping of monopoles starts. When the monopole current is pumped steadily, monopole 
accumulation grows at the edges of the system, inducing a diffusive current. The steady 
monopole distribution is then determined by the balance of this backward diffusion and the 
pumped monopole current. 

Direct evidence of the spin damping monopole is given by detecting the electric current 
discussed above. Surprisingly, the signal from the spin damping monopole might have already 
been detected. In fact, the electric voltage due to magnetization precession has been observed 
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in a junction of a ferromagnet on a Pt film in a pioneering work by Saitoh et al.. 24 The 
mechanism of voltage generation has been explained by the inverse spin Hall effect. According 
to the inverse spin Hall scenario, magnetization precession generates a spin current via the spin 
pumping effect, 25 and the spin current j s is converted into an electric current by the spin-orbit 
interaction (the inverse of the spin Hall effect). This explanation assumes that the conversion 
mechanism of jj = eijkj^j where k is the index representing the spin polarization of the spin 
current. 26 A recent theoretical study has revealed, however, that the conversion formula is 
not universal; it does not apply to the slowly varying magnetization configuration or in the 
presence of disorder. 11 Rather, the conversion formula is an approximated one connecting a 
physical observable (electric current) to a spin current whose definition depends on the specific 
system considered. Our result obtained in the present paper suggests a different scenario that 
is universal owing to the symmetry of Maxwell's equations. 

For the experimental confirmation of the spin damping monopole, of crucial importance is 
the separation of the monopole signal from the inverse spin Hall signal driven by a spin current. 
This is accomplished by applying an electric field (E s ) perpendicular to the junction of Fig. 1. 
The monopole contribution then leads to a transverse electric current as a result of the Hall 
effect of the monopole, 27 while the contribution of the spin current is not affected. In another 
experiment, the strongest evidence of the spin damping monopole is given by observing a 
magnetic field (B s ) produced by monopoles via Gauss's law by electron holography. 

We have shown analytically that a magnetic monopole is a common object in dynamic 
magnetic systems with damping. The control of spin damping monopoles is as feasible as that 
of electrons, for both are governed symmetrically by Maxwell's equations of electromagnetism. 
We here propose the monopolotronics, i.e., the control of monopoles, as a novel concept of 
realizing spintronic devices. 
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